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I . Introduction 

This report covers the period August 1, 1978 to January 31, 1979. 

The primary task during this reporting period was the development of a fully 
conservative solution method and the Incorporation of It Into the TRANDES 
code. Wave drag and massive separation studies were also conducted, 
i I . Discussi on o^f R esearch 

A. Massive Separation Studies 

As Indicated In the last progress report, this effort Is studying the 
application of the present model and method to the NASA 4412 Airfoil. Extensi 
experimental data. Including pressure distributions, have been collected from 
the literature. Unfortunately, difficulty has been encountered in matching 
this data at low angles of attack to theoretical results obtained from tne 
current code. Part of the problem was errors In the coordinate computation 
scheme, but this problem has been corrected and a 1/2° angle of attack error 
still exists. 

Comparison of the fieoretlcal NASA 4412 ordinates with those used 
experimentally Indicates some differences on the lower surface near the 
trailing edge. These discrepancies are small, however, and theoretical 
results using both sets only differ slightly. Thus the origin of the l/2°a 
error Is unknown; and it will for the time being be ignored In the interest 
of obtaining high angle of attack results. These latter results should be 
available shortly, will concentrate on determining input parameters for 
conventional airfoils, and will examine methods of drag calculation. This 
effort is being conducted by an undergraduate student. 

B • Leading E dge & Wav e Drag c^i es 

During this reporting period, work was initiated by an undergraduate 
to develop a leading edge grid Imbedment scheme. Unfortunately, the student 
involved decided in December to accept a job offer and left abruptly. While 
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progress was made, the exact status of the work is currently being studied 
and will, hopefully, be reported later. 

C • ix Conservativ e Codes 

The primary task during this reporting period has been the development 
and Incorporation Into TRANDES of a fully conservative analysis method utilizing 
the artificial compressibility approach. The present v»ork closely follows 
Ref. 1 but allows for lifting cases and finite thickness airfoils and utilizes a 
stretched coordinate system. The solution scheme is a three level Richardson 
method for which has, in the present case, a stability requirement of 




^ r < 


Ml. 


where 


O <. iVCZ 


a "relaxation" like parameter 
ratio of step sizes 

maximum value of horizontal coordinate stretching 
factor 

maximum value of vertical coordinate stretching 
factor 


Initial results with the method Indicated that the solutions trended to 
exhibit ski-ramp shocks. That Is, the pressure just upstream of the shock 
wave would have large changes in dCp/dx. This type of structure is shown on 
Figures 1 and 2 for NVIS *= 1. Numerical experiments at TAMU and NASA Langley 
subsequently showed that this ramp behavior could be mitigated by increasing 
the artificial viscosity just upstream of the shock. In practice, this 
increase was accomplished by using for the artificial viscosity 

yiA = J I 

By increasing NVIS,^would be Increased. 

Typical coarse and medium grid results using the NVIS concept are shown 


on Figures 1 and 2. Obviously, the solution is sensitive to NVIS and, as 
shown by Figure 2, the use of NVIS does not completely eliminate the ski-ramp 
structure. It should be noted that values of NVIS higher than two induced 
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for the medium grid cases numerical Instability. 

At this point it was suggested that the Inclusion of In the 

solution scheme might be of value. For the three-level scheme, the addition 
of means that the governing equation Is of the form 

f 'Pit f ^ f 'Pt f ^ Pxt ' ^ 

where (/)^^ /, ( <P;^ ~ ( ) 

Since the present method utilizes complete surface boundary coiditions and 
hence thickness, the addition of created some problems as to storing 
values at boundary condition ghost points. The creation of a dummy array 
solved this problem, however; and the ensuing results did exhibit enhanced 
numerical stability. 

Simultaneously, NASA Langley discovered that in the ACM formulation being 
used by TAMU that the artificial viscosity was being computed at the (if 1/2, j) 
point Instead of the (i,j) location. A simple fix was devised and typical 
results are shown on Figure 3. Surprisingly, the results were still sensitive 
to the value of NVIS and exhibited for some cases ski-ramp type shocks. 

In spite of these subtle points, the method appeared to be working well. 
Thus, accuracy tests for various cases were conducted. For subcritical lifting 
and nonlifting cases, comparison was made with results from the TRANDFS program. 
For biconvex, the new method, called C0NV3, yielded Cp 

results that agreed almost exactly with the TRAMDES values. 

For lifting cases, such as the aerodynamic coefficient 

agreement was excellent; but the TRANDES results showed Cp overshoot in the 
upper surface leading edge region which was not predicted by C0NV3. Subsequent 
investigation determined that this difference was due to the treatment of the 
(ILE-1 , JB-1) point (see sketch). 




r 


-4- 


The TRANDES code, which solves the full potential equation, needs to compute 
♦ at (lLE-1, JB-1) a 4 value at (ILE,i)B). This value Is obtained by e1th»jr 
satisfying the boundary condition at the lower surface (marked +) or at 
the upper surface (marked X). In TRANDES, the lower extroplatlon was selected 
since this approach yielded better answers, when compared to other methods, 
for the NACA 0012. 

When the TRANDES program was run using the upper surface extropolatlon, 
the results agreed very well with C0NV3. Based upon this agreement. It was 
concluded that C0NV3 was accurate for subcritical cases, although subtle 
leading edge problems might still exist. 

For supercritical cases, comparison was made with non-lifting results 
provided by NASA Langley. Initial comparisons are shown on Figure 4 for a 
10% biconvex case. This case Is difficult In that It has a trailing edge 
supersonic/supersonic shock; and, as can be seen, there Is some disagreement 
between the two methods. Howevfer, the Langley results used a thin wing 
small perturbation boundary condition; while TAMIJ used the full boundary 
condition. When the C?)NV3 code used thin airfoil small perturbation 
boundary conditions, the results shown on Figure 5 were obtained. Subsequent 
tests showed that the differences on Figure 4 were primarily due to usage of 
the ful 1 condition and that finite thickness only influenced the pressures 
slightl’’. In any event, it was concluded that C0NV3 an accurate code, and 
some lifting solutions were obtained. 

At this time (mid December) results of research at NASA Langley by Jerry 
South, Jr. became available which determined the origin of the ski-ramp 
shock problem. As a result a new code, called MIDSE6, has been written. The 
difference, between C0NV3 and MIDSEG will be reported later. Kowever, some 
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results, comparinq the two codes are shown on Figures 6 and 7. Notice that 
the C0NV3 results exhibit a slight skl-ratnp sho*k and that the Cp distribution 
and shock location Is grid size sensitive. On the other hand, MIDSEG does 
not have these features . As a result, the present effort is concentrating 
on the MIDSEG type approach. 

III. Personnel 

(1) Leland A. Carlson, Associate Professor, Principal Investigator, 

1/4 Time August, 1/8 Time September to present. 

(2) Paul Hagseth, Undergraduate Research Assistant, August 1978 
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Figure 3. The Effect of NVIS With Correct/i Location 
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Figure 4. The Effect of Full vs. Small Perturbation Boundary Conditions 
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5. Comparison Of C0NV3 Code Results with NASA Langley 
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Figure 6. Variation With Grid Size 
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